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Line shape analysis of the photoemission intensity from the surface states of Bi2Se3 reveals two 
unusual features: spectral line asymmetry and anomalously enhanced photoemission from the Dirac 
point. The former can be described by the one-particle spectral function assuming that the self- 
energy has an energy-momentum dependent contribution to the imaginary part, which changes 
sign when crossing both the dispersion curves and the energy of the Dirac point. The anomalous 
enhancement can hardly be understood as a self-energy effect, while a final state interference seems 
to be a more plausible explanation. 

PACS numbers: 



I. INTRODUCTION 

The topological insulators, as a new kind of mat- 
ter, and emerging new physics have been predicted 
theoretically^ and later observed in a number of the 
state-of-the-art experiments. Angle resolved photoemis- 
sion spectroscopy (ARPES), applied to the topological 
insulators^ has been playing a major role in their pro- 
motion and study^^. The experimentally observed cone- 
like dispersion of the surface states with the spin de- 
generate Dirac point is now considered as a hallmark 
of the topological surface states^. In this sense, the 
compound Bi2Se3 with a single Dirac cone in the Bril- 
louin zone^'^ provides the most elementary case for the 
in-depth ARPES study of those states. Analyzing the 
spectral weight distribution in spectra of Bi2Se3, we ob- 
serve deviations from the expected one particle spectral 
function: an anomalously enhanced photoemission from 
the vicinity of the Dirac point and unusual asymmetries 
of the surface states spectral function. Having excluded 
the possible experimental artifacts and trivial explana- 
tions, we show that the /ccj-dependent component to the 
self-energy changes sign while crossing both the disper- 
sion curves and the energy of the Dirac point. We also 
conclude that the unexpected intensity enhancement at 
the Dirac point is hard to describe by the single particle 
spectral function. Possible mechanisms for the enhance- 
ments are discussed. 

The paper is organized as follows. In Section II we 
describe the experiment and the results, pointing out 
on the anomalously enhanced photoemission from the 
Dirac point that can be seen in the raw data. In Section 
III, performing a deeper data analysis, we: (i) explain 
why the observed enhancement is unexpected based on 
one particle spectral function and the related sum rules 
(Section III A); (ii) exclude possible artifacts (IIIC); (ii) 
present detailed comparison of the experimental data and 
simple models (HID); (iv) detect an unusual spectral line 
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FIG. 1: The geometry of our ARPES experiment: LH and 
LV stay for the linear horizontal and linear vertical polariza- 
tions of the incident light. 



asymmetry and show that it can be described by one- 
particle spectral function with a certain dependence of 
the self-energy on energy and momentum (HIE), and (v) 
demonstrate that the anomalous enhancement cannot be 
understood if a physically viable model for the self en- 
ergy is considered (III E) . An alternative explanation for 
it, a final state interference, is discussed in Section IIIG. 



II. EXPERIMENT AND RESULTS 

Single crystals of Bi2Se3 were grown using the Bridg- 
man method starting with high-purity Bi and Se ele- 
ments. ARPES experiments have been performed at the 
"1^" and "1^" beamlines equipped with SES 4000 and 
SES 8000 analyzers at the synchrotron radiation source 
BESSY II. Fig. 1 shows the geometry of our ARPES ex- 
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FIG. 2: ARPES image of the Dirac cone in Bi2Se3 measured close to the FK direction, T = 20 K (a). Energy- (b) and 
momentum- (c) integrated ARPES spectra represent partial n(k) and partial DOS, respectively. The difference between 20 eV 
experimental data and the model, shown in panel (b), refers to the right axis. The arrow in panel (c) marks the position of the 
Dirac point. 



periment, in which the entrance slit of the electron an- 
alyzer is fixed, but the sample can be translationally 
moved and rotated. More details can be found in Ref. 7. 
The data have been collected at different temperatures 
in the range from 1 K to 30 K revealing independence of 
the observed anomalies on the temperature. 

In Fig. 2a, we show the aforementioned "hallmark" of 
the topological insulators — a typical ARPES spectrum of 
Bi2Se3 measured close to the FK direction with the sur- 
face state dispersion converging at the Dirac point in the 
center of the Brillouin zone. In contrast to the three- 
dimensional bulk states, which fill here the parabola- 
enclosed region at the top of the spectrum, the surface 
states are truly two-dimensional, which makes them a 
perfect object for an accurate ARPES study of the pa- 
rameters the spectral function depends on. However, we 
have found that the deviations of the photoemission in- 
tensity from the expected spectral function of 2D elec- 
tronic excitations are essential and require thorough in- 
vestigation. 

The most striking effect appears as an anomalously 
enhanced photoemission from the Dirac point. This 
anomaly is clearly seen in Fig. 2b as a small sharp peak 
with the full width at half maximum (FWHM) about 
0.012 A-i on a partial momentum distribution function, 
^Aa;(^), that is defined by integration of the correspond- 
ing ARPES spectrum in the energy range from —0.4 eV 
to —0.2 eV. The range is chosen to be narrow enough 
to get rid of the influence of the bulk states, but wide 
enough to expect a plateau around the center of the Bril- 
louin zone (±0.04 A~^). The local enhancement can also 
be seen (Fig. 2c) on partial density of states, DOSA/e(^), 
that is defined by integration in the momentum window 
±0.12 A~^. In order to understand why the enhance- 
ment is anomalous, one needs to consider the structure 
of the spectral function and its relation to the ARPES 
spectrum. 



III. DATA ANALYSIS AND DISCUSSION 

A. Trivial spectral function 

The photoemission intensity from 2D or quasi-2D 
states are expected to be well described by the product 
of the Fermi function, /(cj), and the one-particle spectral 
function, A(co',k) (a function of energy and in-plane mo- 
mentum). For the photoemission current from the n-th 
band, neglecting the experimental resolution and photoe- 
mission background, one can write 

/n(c^,k) = M(n,/ii/,cj,k)/(cj)An(a;,k), (1) 

where the transition probability (the squared modulus of 
photoemission matrix elements) and some extrinsic fac- 
tors (such as detector efficiency) are included in a fac- 
tor M, and An is supposed to have a simple represen- 
tation in terms of bare electron dispersion, e'^ of the n- 
th band, and the corresponding one-particle self-energy. 

For the sake of simplicity, we will consider the spec- 
trum such as shown in Fig. 2a, i.e. measured in a cut 
of the k-space through the center of the Brillouin zone, 
as one-dimensional, i.e. k = {kx^^) = {k^O). In Ap- 
pendix A we show that the finite momentum resolution 
along ky does not affect the eligibility of this approxima- 
tion. Then, the spectrum can be considered as a sum 
of intensities from two bands, hi^l) + ^2(^^)7 the states 
in which can be enumerated in a way to trace smooth 
monotonic dispersions ej. = = Ek- So, one can write 

I{uj,k)=h{uj,k)^h{uj,-k). (3) 
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This means that one does not expect to see any sharp fea- 
tures in energy or momentum integrated /(cj, k) except 
those present in similarly integrated Ii{uj^k). 

The dependence of the coefficient M on k and uj in 
Eq. (1) is determined by the photoemission matrix ele- 
ments and can be neglected in our analysis. Indeed, the 
enhancement is localized within less than 0.05 eV, which 
is by an order of magnitude smaller than one would ex- 
pect from the variation of matrix elements with energy^. 
In addition, we have observed the effect at several exci- 
tation energies, as shown in Fig. 2b, c. Both observations 
practically rule out M(co', k) as a source of the anomaly 
reported here and, in the following, we focus on the one- 
particle spectral function and examine whether the ob- 
served enhancement can be described by a sophisticated 
cjk-dependent self-energy or should be caused by an ex- 
trinsic mechanism. 



B. ARPES sum rules 

Now, let us remind the sum rules that the one-particle 
spectral function must comply with.^ The simplest one 
comes from its definition and says that its integral over 
all frequencies is unity: 

/CO 
A{uj,k)duj = 1. (4) 
-oo 

A more practical one, 

/oo 
A{uj,k)f{uj)duj = n(k), (5) 
-oo 

gives the momentum distribution function n(k) that is 
close to 1 when A{uj) is peaked well below the Fermi- 
level, i.e. when (^k + ^')^ > ^"'^ for + < 0. 

A similar function, which can be extracted from 
ARPES spectra, we will call a "partial" momentum dis- 
tribution function, 

nUJ2 

nAc.(k)= / A{uj,k)f{uj)duj, (6) 

where Auo = 002—001 is the range of integration that is 
always finite in experiment, but in this case we narrow it 
deliberately to avoid the influence of the electronic bands 
that are out of interest. 

The difference between n(k) and nAu;(k) is due to the 
tails of A{uo) that are cut by the window of integration. 
Since, within the assumptions made above, the profile 
shown in Fig. 2b represents uau (k) , one can expect that 
the observed peak might be associated with a sharp nar- 
rowing of the energy distribution curves (EDC) at /c = 
that, in turn, corresponds to a sharp minimum in S''(k) 
at this point. 

Similarly, while the summation over all k-states gives 
the density of occupied electronic states 

DOSH = ^A(w,k)/H, (7) 

k 



it is often useful -^^'-^-^ to define a partial density of states 

BOSAkH = J A{uj,k)f{uj)dk, (8) 

which can be compared to the experimentally derived 
functions such as shown in Fig. 2b. In our case, the re- 
gion of integration is determined by the momentum win- 
dow Ak from -0.12 A"^ to 0.12 A"^ along the slit of the 
analyzer and by the momentum resolution perpendicular 
to the slit, the effect of which is considered in Appendix 
A. 

One can see from Eq. (2) that both DOS(cj) and 
DOSa/c(^) are roughly constant when Sk is close to linear 
and S can be considered as weakly k-dependent. This 
has been shown to be the case for a number of mate- 
rials including strongly correlated cuprates.^^'^^ Thus, 
again, the observed enhancement of the photoemission 
signal seen in Fig. 2c at the Dirac point energy might be 
caused by a sharp narrowing of the momentum distribu- 
tion curves (MDC), which is a consequence of a sharp 
minimum in Ti^\oo) at this energy. 

In Section HIE we examine whether it is possible to 
describe simultaneously the peaks on both nAu;(k) and 
DOSa/c(^) functions with a cjk-dependent self-energy, 
but before, we have to prove that the effect is real, ex- 
cluding possible experimental artifacts. 



C. Excluding trivial artifacts 

First, one may question whether such an effect can be 
a result of an inaccurate position in k-space or of a finite 
momentum resolution. Both cases are similar to the case 
of a finite gap due to lifting of the degeneracy at the 
Dirac point. Any of these reasons would result in a dip 
rather than a peak in DOSAfe('^) and should not affect 
n/\uj{k) at all (see Appendix A) which is in contrast to 
our observation. 

Second, one should examine a possible non-linearity of 
the detector, which has earlier been reported to be an is- 
sue with similar analyzers^^^^. Indeed, the observed en- 
hancement can be qualitatively understood if the readout 
intensity as a function of actual photocurrent increases 
faster than linearly. One can notice, however, that it 
is an unlikely reason for the observed anomaly since its 
width in momentum is very narrow, less than a half of 
the minimum line width that can be found in the spec- 
trum. Nevertheless, in Appendix B, by careful analysis 
and modeling the ARPES spectrum, we show that the 
actual non-linearity of the detecting system we use can- 
not result in the observed anomaly. 



D. Fit to the model 

In order to better reveal the differences between the 
measured spectrum and the expected model further, we 
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FIG. 3: The experimental data (a), the result of the fit-based model (model 1) (b), and their difference (c), superimposed by 
the MDC dispersion (dashed lines). The spectral function parameters derived from the MDC fit, as described in the text (d): 
Ti" is an effective imaginary part of the self-energy, M — vM^j, where Muj is the MDC area and v is the bare band velocity, 
B is the momentum independent background in counts per channel. Examples of the MDC fit at 0.275 eV (e) and 0.302 eV 
(f) binding energies: experimental data (black symbols), two-lorentzian fit (red solid line), and their difference (dashed blue 
line). The experimentally derived band dispersions and the band velocity (g). Partial momentum distribution nAa;(k) (h) and 
partial DOS (i) derived from experiment (solid black line) and the model 1 and 2 as well as for the case of 10% nonlinearity 
applied to the model 1. 



consider the spectrum as a set of MDC's^^ and fit them 
to 

/(cj, k) = h{uj, k) + /i(a;, -k) + B{uj), (9) 
rewriting /i(co', k) in an MDC form: 



Ii{uj,k) 



W{uj) 



[km{0j)-k]^^W^{uj)' 



(10) 



Here = M/v; the MDC width (half width at half 
maximum) W{uj) = —T,"{uo)/v{uo)] the position of MDC 
maximum kmioj) = [uo — £q{uo) — Y^' {ijo)]/v{uj)] £q = 
ek{]^m) - v{uj)km; v{uj) ^ (dkm/duj)-^; B{uj) stays for 
both the momentum independent background and con- 
tribution from the bulk states. 

In Fig. 3 we shows the results of such a fit and compare 
the experimental spectrum (a) to the model (b). Panel 
(c) represents their difference. The two Lorentzian fit 
to the experimental MDC is stable everywhere except 



± 5meV around the Dirac point. So, in the first step we 
derive km{^) in the fit stability range, interpolate it in the 
vicinity to the Dirac point and smooth it with the 10 meV 
smoothing window. The result is show in Fig. 3g. Then, 
holding km{^)^ we derive other parameters, using the 
only B > {) constrain, and put in the model (model 1) the 
smoothed functions T^"{uo)^ M{uj)^ and B{uj)^ as shown 
in Fig. 3d. By the smoothing, we filter out the rapid 
variation of the fitting parameters in the close vicinity to 
the Dirac point. 

A one-to-one comparison of the experimental and 
model MDCs shows that the experimental spectrum dif- 
fers from the expected one in two ways. First, an unex- 
pected MDC lineshape is observed in the vicinity to the 
Dirac point: a narrow peak has appeared in a middle of 
two nearly merged peaks, as shown in Fig. 3f. Evidently, 
it is this peculiarity that is responsible for the afore- 
mentioned enhancement of the photoemission intensity 
in the energy- and momentum-integrated spectra. Sec- 
ond, there are systematic deviations of the experimen- 
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FIG. 4: A pattern of the energy dependent part of the self- 
energy: (a) sign of S{(jU, k); (b) sign of a" {u^ k), (c) suggested 
symmetry (center), and (d) their difference; positive is red, 
negative is blue, zero is white. 



tal MDCs from the line shape expected for the spectral 
function with the momentum independent self-energy, as 
shown in Fig. 3e. Those deviations can be understood 
as an asymmetry of each MDC peak in respect to the 
position of its maximum. Fig. 3c shows that these devia- 
tions are systematic. In the following sections we discuss 
the both anomalies separately, but first, having in hands 
the developed fit-based model, we address the question 
of detector linearity more carefully. 

Assuming a weak non-linearity of the the read-out 
count rate on real photoemission intensity D = /(I + 
al /Imax)i we can show that for the detector and count 
rates we use, a < (see Appendix B), i.e. one would 
expect a depletion rather than enhancement of the in- 
tensity in the Dirac point. Nevertheless, we model the 
case of a = 0.1, recently reported for similar detectors^ ^. 
In Fig. 3(h,i), one can see that although the effect of 10% 
non-linearity is noticeable, it cannot produce such sharp 
peaks as observed in nAoj{k) and DOSa/c('^)- 



E. Line shape asymmetry 

Let us assume that the difference between the real and 
model spectral functions, such as shown in Fig. 3c, is 
the result of a small component of the self-energy with 
strong dependence on energy and momentum, and find 
its simplest form that can explain both the spectral line 
asymmetry and the anomalously enhanced photoemis- 
sion from the Dirac point. 

We start with the asymmetry issue and write the total 
self-energy as S(co',k) = So[l + cr(cj,k)], where is 
a main self-energy part that may be constant or may 



vary slowly with uj. For <C 1, completely neglecting a 
contribution of to the spectral function (that violates 
the Krames-Kronig consistency^^ but does not affect our 
finite result), one can derive: 



AA{uj, k) = A(S) - A(So) - Sa", 



where 



^0 



^2 



^0 



y" 



(11) 



(12) 



and = w - (Tfc - 
i.e., for l^fc - 
e.g., for a < E'o, 



- Eg. Far away from the Dirac point, 
» 2Sq, and close to the dispersion. 



AA{co, k) = a"{co, fc)/7rS^'. 



(13) 



From Eq. (11) one can see that AA{uj^ k) changes sign 
when either S{uo^k) or (j"{(jj^k) does. The former hap- 
pens when the spectral function intensity is about half 
of its maximum, see Fig. 4a or Eq. (12), but the sign of 
AA(a;, k) changes also at ^/^ = 0, = 0, and uj = co^df, 
as it follows from Fig. 3c or Fig. 4b. One can conclude 
therefore that a" should change the sign when crossing 
either the dispersion curves or the energy of the Dirac 
point, as it is shown in Fig. 4c, which can be a signature 
of a spin- flip scattering^^. 

Note that while the sign-change pattern of a" shown in 
Fig. 4c can be described by its dependence on uj only, the 
dependence on k is still necessary to explain the asym- 
metry of the MDC peaks. 



F. Anomalous peak at Dirac point 

One can certainly describe a single MDC with a mo- 
mentum dependent self-energy. The sharp peak that ap- 
pears in the middle of the MDCs in the vicinity to the 
Dirac point, as shown in Fig. 3f, would require the scat- 
tering rate to have a sharp drop at /c = 0. In Fig. 5 (top), 
we fit the model MDC, defined by Eq. (9), to the experi- 
mental one writing the imaginary part of the self-energy 
as Y."(k) = J:'^[l^S'\k)], where S"{k) = awl/iP^wl), 
with = -0.024 eV, a = -0.4, Wk = 0.0011 A-\ 
km = 0.0068 A~^, and the band velocity at the Dirac 
point V = 2.12 eVA. 

However, further analysis shows that these parame- 
ters do not describe the whole spectrum and, more- 
over, that the dependence of the self-energy on k only 
cannot describe both n/\^{k) and DOSa/c('^) simultane- 
ously. In particular, fitting the model to experimental 
nAuj{k) gives very different Sq = —0.02 eV, a = —0.7, 
Wk = 0.008, and, in addition, results in three times higher 
peak in the modeled DOSAfe(^) than it is observed ex- 
perimentally. 

One may hope that the problem can be resolved if S" 
depends on both k and cj, but we have found that it can- 
not. Despite the linear dispersion, the influence of the 
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momentum and energy dependences of the self-energy 
on nAuj{k) and DOSA/e(^) is not symmetric. The reason 
for this is in the Kramers-Kronig relation between en- 
ergy dependences of the real and imaginary parts of the 
self-energy. One can show that if the Kramers-Kronig 
consistency is taken into account, the energy dependence 
of Ti" has virtually no effect on n/\^{k). To illustrate 
this, we model the depletion of T^" as 



a 



1 



(14) 



This assumes the Kramers-Kronig consistent counterpart 



5\uj,k) = 



(15) 



contributes to the real part of the self-energy as 
T.'{uj,k) = T.'i5'{uj,k). 

In Fig. 5 (bottom) we show that both n/\uj{k) and 
DOSa/c(^) can be fitted simultaneously only in an ar- 
tificial case when S' is set to zero. This case is pre- 
sented by the solid curves which are modeled to fit to the 
corresponding experimental n/\^{k) and DOSa/c(^) with 
the following parameters: Hg = —0.02 eV, a = —0.3, 
= 0.04 eV, Wk = 0.008 A~^. One can see (dashed 
curves) that taking into account the Kramers-Kronig 
consistency essentially eliminates the peak in n/\^{k). So, 
the energy dependence of T^" cannot help to reconcile the 
discussed anomaly in the experimental spectrum with the 
one particle spectral function. This calls for more rigor- 
ous consideration of the photoemission process. 



G. Final state interference 

A plausible explanation for the observed enhancement 
and interference-like spectrum might be a finite state in- 
terference of the photoelectrons emitted from the k^^ and 
— states in the vicinity to the Dirac point. 

Let us assume that (1) the width of a single MDC is 
defined by the size, (i, of the region from where the pho- 
toelectrons are emitted coherently, and (2) the photoe- 
mission from k^ and — k^ states from this region remains 
also coherent. For the sake of simplicity, the former can 
be described by the Fraunhofer diffraction and for the 
photoemission intensity from a single /c^ state one can 
write 



/i(a, 0) (X 



sin[^ sin(a + 0)] 



^ sin(a + 0) 



(16) 



where a is the angle of electron emission, is the angle 
which corresponds to the initial state of electron, such as 
hkuj = ^/2mEk sin(^), ^ = iid/X^ is the kinetic energy 
of photoelectron, A = 27:h/ ^2mEk is its wave length, 
and m is the electron mass. 

Then, based on the second assumption and for small 
a and ^, the photocurrent from k^ and —k^ states can 
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FIG. 5: Momentum dependent self-energy. The peaked 
MDC can be a consequence of a sharp decrease of the scatter- 
ing rate at /c = (top). Both nAuj{k) and DOSAfc(<^) can only 
be fitted to the experiment if the Kramers-Kronig consistency 
is not taken into account (bottom). 



be written, instead of Eq. (3), as 



/(a, 6) (X 



sin[^(a + 6>)] sin[^(a-6>)] 



(17) 



that is similar to the famous two-slit diffraction pattern. 

In Fig. 6 we show both, /(a, 6) for 6 = 0.02 rad, and a 
function 



n{a) 



I{a,0)dO 



(18) 



that is analogous to the momentum distribution function 
discussed above. Ej^ = 15.5 eV and ^ = IOOtt that corre- 
sponds to A = 3 A and d = 300 A. The parameter ^ has 
been chosen to reproduce the width of the experimental 
MDCs. 

As one can see, the simple interference model can ex- 
plain both the appearance of the k = peak in n{k) 
and the fact that it looks narrower than the MDC peaks. 
So, although the photoemission from k^j and —k^j states 
is usually considered as non-coherent, one may assume 
that a finite electron-electron interaction can set the re- 
quired coherence while the extended nature of topolog- 
ical surface states^^, due to very long phase-breaking 
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lengths^^, may provide a sufficient area on the sample 
surface from which the coherent photoemission occurs 
and results in the observed effect. The coherence of the 
and —k^ states can be a consequence of a coherent 
spin rotation, recently suggested to explain the interfer- 
ence of spin states in photoemission from the surface al- 
loy Sb/Ag(lll)^^. Among other examples of interfer- 
ence in photoemission, that looks similar to our case, we 
want to mention the coherent photoemission from dimer 
molecule^^, or the spin-state interference, observed in res- 
onant photoemission induced by circularly polarized light 
from magnetized Gd^^. 



them to the case of a finite band gap A due to lifting of 
the degeneracy at the Dirac point. We model the data 
as 



/2 



(Al) 



dh 



where Rk is the momentum resolution perpendicular to 
the slit of analyzer, i.e. along ky^ 



A 

2"' 



£± = ±Vf Jkl + k^- WDP, 



(A2) 
(A3) 



is the bare band velocity, cjdp is the energy of the 
Dirac point. Here we use = 2 eVA, 0.32 eV, 

and T^" = 0.02 eV, to be close to the experimental data. 

The images in Fig. 7 represent the examples of 
^mod(^7^x) calculated for different values of each of the 
parameters A, ky^ or Rk^ while others where kept zero. 
One can see all these cases result in a dip rather than 
a peak in DOSAfe(^) and, by definition, do not effect 
^Au{k)' This result is natural since the dispersion at fi- 
nite ky is similar to the dispersion with finite A, which 
is known (e.g., from tunneling spectroscopy) to result in 
a dip in DOS. Consequently, at finite momentum resolu- 
tion, one cannot expect different shape of DOS, which is 
just an average of partial DOSs with different ky. 



IV. CONCLUSIONS 

We have found the profound deviations of ARPES 
spectrum of the topological surface states of Bi2Se3 
from a generally expected spectral weight distribution 
of weakly interacting quasiparticles. The observed spec- 
tral line asymmetry can still be described by the one- 
particle spectral function but with the self-energy that 
strongly depends on both energy and momentum. The 
kcj-dependent constituent of the imaginary part of the 
self-energy that changes sign when crossing both the dis- 
persion curves and the energy of the Dirac point can be 
a consequence a spin- flip scattering^^. In opposite, the 
the anomalously enhanced photoemission from the Dirac 
point is unlikely to be a self-energy effect if the Kramers- 
Kronig consistency is taken into account. In this situa- 
tion a final state interference looks more plausible expla- 
nation for it. 



Appendix A: Effect of finite band gap and 
momentum resolution 

In Fig. 7 we examine the effects of not exact position in 
k-space and a finite momentum resolution, comparing 



Appendix B: Detector linearity 

In order to estimate the linearity of the detector we com- 
pare pairs of standard ARPES images, such as shown 
in Fig. 1, measured with different photon flux. Let us 
express the read-out intensity I) as a function of real 
intensity / in each detector channel as D^'^ = D{Iij). 
Approximating the detector function as !)(/) ~ / + a/^, 
and assuming that real photoemission intensity can be 
changed from / to nl by changing the photon flux with- 
out changing other parameters of the light (such as en- 
ergy, polarization, etc.), for two spectra measured with 
different photon flux one can write 

Di = D{I)^I{l^aI), (Bl) 
Dn = D{nl) ^ nl{l + ani). (B2) 

Then, for small a/, 

Dn ^ nDi + an{n - 1)D^, (B3) 

so, one can estimate both n and a comparing two spectra 
measured with different photon flux. 

In Fig. 8, we plot D^i^ vs DW i.e., the intensity of 
each pixel of one ARPES spectrum vs the intensity 
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FIG. 7: Simulation of a finite band gap (top row), not exact position in k-space (middle row), and a finite momentum 
resolution (bottom row). All images are shown in the same absolute color scale. The dispersion, if shown, is presented by a 
white dotted line over corresponding images. 



of the corresponding pixel of another spectrum. Here 
we use the data for three spectra: #1, #2, and #3, 
that differ by size of the photon beam from the un- 
dulator (controlled by the undulator blades: 0.5x0.5, 
0.8x0.8, and 4x4 mm^, respectively). So, we plot #2 
vs #1 data in Fig. 8 (left panel) and #3 vs #2 data in 
Fig. 8 (right panel). The data are shown in counts per 
pixel that is 1.5meVx0.02° of size and is one of 294192 
(454x648) pixels in each image. For the recorded spec- 
tra the time per spectral channel was 453 s, so the max- 
imum counts recorded per pixel, Imax = 6-10^ cpp, cor- 
respond to the count rate 13 counts per pixel per sec- 
ond, or to 0.43 counts per second per 1 /im^ of the de- 
tector area. The fit of the data in Fig. 8 (left panel) 
to Eq. (B3) gives n = 2.509 ± 0.002 and negative a = 
— 1.04 X 10~^ ± 1.5 X 10~^ cpp~^, or, in a dimensionless 
form, a = almax = —0.165 ± 0.002. For comparison, we 
add three lines that represent Dn{Di) function for dif- 
ferent a = 0, 0.1, and 0.4. Those values are chosen to 
be close to the experimentally measured non-linearity of 
different detectors: 10% for SPECS Phoibos 150^^ and 
40% for the Scienta SES200 detectors^^'^"^. The fit of 
the data for higher count rates (Fig. 8 right panel) gives 
n — 6.5 and a = —0.05. 

Since such a procedure of the evaluation of the de- 



tector linearity is not widely used, we demonstrate its 
applicability in Fig. 8 (central panel). We model both 
Di and Dn based on the spectrum #1, taking the ex- 
perimental data as / and calculating D\'^ = D{Iij) and 
Dl{^ = D{nli-^ij-^i) with n = 2.5 and different a = 0, 
0.1, and 0.4. The one-pixel shift of Dn is introduced to 
simulate the experimental noise. The fits of the corre- 
sponding distributions to Eq. (B3), shown for each dis- 
tribution, are also shown in the left panel for comparison 
to the real data. 

This analysis allows us to conclude that the linearity of 
the new Scienta detectors is significantly improved: the 
non-linearity of the detector we use is of the opposite sign 
but much smaller in absolute value than that has been 
reported for the old detectors^^'^^. 
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FIG. 8: The intensity of each pixel of one ARPES spectrum, Dn{(^i, kj), vs the intensity of the corresponding pixel of another, 
Di, spectrum: spectra #2 vs #1 (left); #3 vs #2 (right); and simulated based on spectrum #1 with n = 2.5 and a = 0, 0.1, 
and 0.4 (center), as described in the text. 
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